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Abstract 

We address how to construct an infinitely cyclic universe model. A major consid¬ 
eration is to make the entropy cyclic which requires the entropy to be re-set to zero 
in each cycle expansion —> turnaround —> contraction —> bounce —> etc. Here we 
re-set entropy at the turnaround by selecting the introverse (visible universe) from 
the extroverse which is generated by the accelerated expansion. In the model, the 
observed homogeneity is explained by the low entropy at the bounce, The observed 
flatness arises from the contraction together with the reduction in size between the 
expanding and contracting universe. The present flatness is predicted to be very 
precise. 


1 Introduction 


It is of broad interest among all physicists, indeed all scientists, to understand the early 
universe especially whether a singular bang ever occurred as suggested by the most simple- 
minded interpretation of the Friedmann equation. The idea of a big bang has proved 
remarkably resilient despite the knowledge that extrapolating back the Friedmann equation 
to the singularity must somehow be flawed. 

The discovery [1] of the cosmic microwave background(CMB) resolved a dichotomy then 
existing in theoretical cosmology between steady-state and big-bang theories. The inter¬ 
pretation of the CMB as a relic of a big bang was compelling and the steady-state theory 
died. Actually in 1965 it was really a trichotomy being reduced to a dichotomy because a 
third theory, a bounce in a cyclic cosmology, had been under study since 1922 |[2J. 

An infinitely cyclic cosmology was always the most popular alternative to an initial singu¬ 
larity but its correct development was hindered by an impossible-seeming stumbling block 
created by the monotonically increasing entropy of the universe. Entropy is sometimes 
perceived as not a fundamental concept because it involves a very large number of particles 
and is therefore not a property of a single particle. Nevertheless, in studying theoretical 
cosmology we shall argue that the entropy of the universe is the single most important 
concept. Interest in the pursuit of a successful cyclic model is reflected in the existence of 
three recent popular books [3H5], as well as in [6]. 

Nevertheless, for purely theoretical reasons, the cyclic model had been discarded due to 
the Tolman Entropy Conundrum(TEC) [7JB]- The TEC, stated too simply, is that the 
entropy of the universe necessarily increases, due to the second law of thermodynamics, 
and therefore cycles become larger and longer in the future, smaller and shorter in the 
past, implying that a big bang must have occurred at a hnite time in the past. 

Some progress towards a solution of the TEC was made in [9] using the Come Back 
Empty (CBE) assumption in the so-called BF model. A huge entropy was there jettisoned 
at turnaround and the significantly smaller universe, empty of matter, contracted adiabat- 
ically to a bounce with zero entropy. The original BF model employed so-called phantom 
dark energy with equation of state u> < — 1. Since this violates every well-known energy 
condition, it appears more sensible to employ a cosmological constant with w = -1 as we 
insist upon here, in a standard A CDM model of the present era . We shall find that CBE 
at turnaround is still possible. An abbreviated presentation of some of the present results 
appeared in in. 

The emphasis in this article is on entropy and the second law of thermodynamics. No 
cause or mechanism will be provided for the turnaround from expansion to contraction of 
the bounce from contraction to expansion. The only theory used will be general relativity, 
in particular the Friedmann equation to describe the expansion and contraction era of 
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the space-time described by a homogeneous and isotropic FRLW metric. At almost all 
times general relativity is sufficient. For the turnaround and bounce, an unknown extra 
ingredient is necessary as a part of a C-theory where C means Complete. Our assumption 
is that use of the Friedmann equation is sufficient to draw conclusions about the cyclic 
cosmology. 


To illustrate a phenomenological fudge factor to precipitate the turnaround and bounce 
one could modify the Friedmann equation to 


a(t) 
a(t ) 


8 v tG 


PTorif) 




k (t) 

a(t ) 2 


( 1 ) 


where ty, ts are the turnaround, bounce times respectively and where the fudge factor in 
square brackets has no derivation and violates diffeomorphism invariance so is merely to 
illustrate a point, in the absence of any known C-theory. In mi it is shown that infinite 
cyclicity requires that tr ~ 1.3 Ty, at which time the radius of the extroverse is Rev^t) ~ 
10 67 m, a length scale far beyond where general relativity has been tested. At the bounce, 
if it is at the Planck time, the radius of the extroverse is Rev^b = 10 - 44 s) ~ 10 -6 m, 
which is a length scale where Newton’s law of gravity has not even been tested. 


Eq. (Q]) is a good approximation to the Friedmann equation 



87 tG 

~1T 


pTOT\t ) 


k (t) 

a[t ) 2 


( 2 ) 


for all times except very close to t ~ or t ~ because a(t T ) is extremely large and 
a(ts) extremely small. In this paper we shall use only Eq.(j2]) and make the reasonable 
assumption that corrections from effects of C-theory near to the turnaround and bounce do 
not significantly change conclusions about entropy and the second law of thermodynamics. 


Most of research on cosmology and the early universe in the last few decades have assumed 
an inflationary paradigm where there is a short burst of superluminal accelerated expansion 
early on in the universe. The present cyclic model does not use any such assumption. An 
inflationary stage is redundant because the geometry becomes flat without it. This is 
because the value total( t) = 1 is a stable fixed point of the contraction equation, just 
as it is an unstable fixed point under expansion. This fact alone is sufficient to achieve the 
observed flatness without inflation. The CBE model makes the flatness much more exact 
to a degree that the departure from IFtot^o) = 1 at the present time is too small to be 
observable. 


The plan of this paper is that in Section 2 we discuss the present expansion era; in Section 
3 we study the turnaround from expansion to contraction; in Section 4 the topic is the 
bounce from contraction to expansion; finally, in Section 5, there is a concluding discussion. 
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2 Present expansion era 


Assuming the cosmological principle of homogeneity and isotropy 
metric 


ds 2 = dt 2 — a{t) 2 


dr 2 


1 — k(t)r‘ 


+ r 2 (d,9 2 + sin 2 9d(f) 2 ) 


leads to the FLRW 

(3) 


where a(t) is the scale factor and k the curvature. Inserting this metric into the Einstein 
equation leads to two Friedmann equations. The first is the expansion equation 



87 tG 


PTor(t) 


k (t) 

a{t) 2 ’ 


(4) 


where pTorif) is the total density. 


Using the continuity equation p(t) + 3 H(t)(p(t) + p(t )) = 0 , differentiation of Eq.fllD gives 
rise to a second equation 


a(t) 

a(t) 


4t tG 

~1T 


(p(t) + 3p(t)). 


(5) 


The critical density is defined by p c {t) = [3H(t) 2 /8 ttG\ and discussion of flatness involves 
the proximity to unity of the quantity 


a 


TOT 


(t) = 


Ptot 

Pc(t) 


( 6 ) 


Let the present time be t = to — 1.38 x 10 10 and normalize a(t) by a(t 0 ) = 1. Other useful 
cosmic times, all measured relative to the would-have-been bang at t = 0, are the Planck 
time at tpi anc k = 10 ~ 44 s, the electroweak phase transition time at tpw = 10 - 10 s, the onset 
of matter domination at t m = 4.7 x 10 4 y and the onset of dark energy domination at 
Oe = 9.8Gy. 

In Table 1 we summarize the well-established history of the present expansion era. 
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Table 1: History of present expansion. 


TIME 

t 

Scale factor 
a(t) 

Comments 

13.8Gy 

1 

present time 

9.8 Gy < t < 13.8 Gy 

a(t) = exp[i7 0 (t — to)] 

dark energy dominated 

top = 9.8 Gy 

0.75 

onset of dark energy domination 

47 ky < t < 9.8 Gy 

a(t) oc t 2/3 

matter dominated 

47 ky 

2.1 x 10 -4 

onset of matter domination 

t < A7ky 

a(t) oc t l/2 

radiation dominated 

t = 10~ lo s 

1.7 x 10~ 15 

clectroweak phase transition 

icr 44 s < t < io~ 1(J s 

a(t) oc t 1/2 

possible inflation or bounce 

O 

1 

Co 

1.7 x 10- 32 

Planck time 


In the unadorned big bang theory, one has dependences of the scale factor a(t) f 1//2 for 
t Planck < t < t m , a(t ) ~ t 2/3 for 4t m < t < t DE , and a(t) = exp [H 0 (t - t 0 )] for t DE < t. 
This leads to the values of a{t) displayed in Table 1. 

The two most striking observations of the visible universe are isotropy to an accuracy 
l±O(10~ 5 ) and [12| flatness tt T oT(t o) = 1.00 ±0.05. In big bang theory the latter implies 
that in the early universe Otot( 1) is given at t — t E w ,tpianck by 

I^Itot^ew) — 1| ~ 0(10 26 ), (7) 

\^TOT(tPlanck) = 1 | ~ 0(10 6 °), ( 8 ) 

as we shall re-derive later by considering the contraction of a cyclic universe. Once we 
incorporate cyclic entropy, the resultant flatness will be even more extreme at the Planck 
time than shown in Eq.(JH]), see e.g. Eq. fl3ip below for an example of this phenomenon. 

We shall focus first on the geometrical condition of flatness, then return towards the end 
to address homogeneity in terms of, what is an intimately related concept, the low entropy 
of the universe at the bounce. 

A convenient way to discuss flatness is to rearrange Eq. (j4j) after division by 77(f) 2 as 

(SW«) - 1) = (9) 


In a decelerating expansion the denominator of the RHS in Eq.Q becomes smaller and 
smaller and Htot(^) deviates more and more from Otot( 1) = 1- This is why the proximity 
of [Htot(^o) — 1] to zero imposes the strong initial conditions in Eqs.([7D and (JSJ). 
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The most popular explanation of flatness is “inflation” P3HEZ! which can be defined as 
insertion of a brief period of highly accelerated superluminal expansion at a time t ~ 
^inflation during the radiation-dominated era t P i anck < t in fiation < t ew■ While inflating, 
a(t) becomes extremely large enforcing flatness sufficiently so precisely in Eq.(l9l) that the 
subsequent decelerating expansion during ti n fi a uon < t < toE does not remove it. 

Inflation also explains homogeneity and has other successes including the prediction of 
scale-invariant density perturbations. It has been claimed that inflation has been success¬ 
fully incorporated into string theory [T8lfl9]. The successful discovery of the BEH boson 
adds credibility to the existence of the scalar inflaton. One possible objection to inflation, 
that it leads to eternal inflation ra and hence to a multiverse in which predictivity is 
hampered by the measure problem GU. is not a fatal flaw. Only a compelling alternative 
theory could cast doubt on the correctness of inflation: cyclic entropy can provide such 
an alternative as discussed in the remainder of this paper. 

If one insisted on criticizing inflation (any theory can be criticized!), it would be that the 
initial conditions before inflation are unspecified. 
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3 Turnaround from expansion to contraction 

To solve the TEC, the entropy of the introverse immediately after turnaround must es¬ 
sentially vanish as discussed in [9]. As an aid to discussion of cyclic entropy, we display 
Tables 2 and 3 which illustrate the future evolution of A CDM models. 


In Table 2 we start from the present time t = f 0 = 13.8 Gy and proceed into the future first 
in steps of 10 Gy, then 50 Gy, then 250 Gy and finally 500G5/. In the second column is the 
corresponding value of the scale factor a(t) with the conventional normalization a(to) = 1. 

In the third column, we assume the number of galaxies inside the introverse at present, 
t — to, is exactly one trillion, 10 12 . Any departure would involve a corresponding rescaling 
of all the entries in this column, but the entries < 1 will be unchanged by any reasonable 
such departure. 

For the fourth column in Table 2, the radius of the introverse (IV) is calculated as usual 
from ^ 

Riv(t) = AAGly + c I dt a(f) _1 , (10) 

Jt 0 

with a(t) = exp [H 0 (t — to)]- The value of Riv(t ) gradually evolves to an asymptotic value 
of Riv{ oo) = 58 Gly, limited by the speed of light. 

Because of the superluminal expansion of space, Rjv(to) is increased to the much larger 
radius Rev(P) = a(t)i?/y(t 0 ) by the stretching of space, as shown in the fifth columnn of 
Table 2. We refer to this larger space as the extroverse (EV). 


It now requires very great care to identify the appropriate introverse at turnaround to 
fulfill [9] the CBE assumption necessary for cyclic entropy. One might initially be tempted 
to include the Milky Way galaxy which we inhabit. That would be too Ptolemaic and a 
most grievous error! The correct choice of introverse at t — tr is instead a sphere which 
contains no matter at all, luminous or dark, including no black holes. It contains instead 
only dark energy with no entropy and small quantities of both curvature and radiation, 
this last being inevitable and actually crucial to the ensuing derivation of flatness. 

In the discussions of [9j, a central role was played by the fraction (/) defined by 

= Riv(t T ) _ Riv{t T ) 

a(tT)REv(t o) Rev^t) 

This all-important quantity is listed in the sixth and last column of Table 2. 


Because the present time t — to is not special, we have taken the liberty to display Table 
3 which is a slightly different but more physical version of Table 2. Table 3 begins at the 
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physically important time t = toE = 9.8 Gy which is the onset of dark energy domination 
and hence of the superluminal accelerated expansion. We highlight the second row which is 
the present time t = t 0 = tpE+^Gy, because it illustrates that aextroverse presently exists 
with Rev = 52Gly and volume ~ 60% larger than the introverse. This must contains some 
600 billion additional galaxies which will remain forever unobservable, given a continued 
accelerated expansion. 


Table 2: Dependence on tr normalized to Rivit o) = Rev it o)- 


Turnaround 
(hr — to) 

a(tr ) 

Galaxies 
in introverse 

Rivitr) 

Rev^t) 

f 

0 

1 

o 

to 

UGly 

UGly 

1 

lOGy 

2.1 

1.8 x 10 11 

51 Gly 

91 Gly 

0.56 

20Gy 

4.3 

2.4 x 10 1U 

55 Gly 

190 Gly 

0.29 

30Gy 

8.8 

3.0 x 10 9 

56 Gly 

390 Gly 

0.14 

40Gy 

18 

3.6 x 10 8 

57 Gly 

800 Gly 

0.071 

50Gy 

37 

4.8 x 10 7 

58 Gly 

1.6 Tly 

0.036 

lOOGy 

1.5 x 10 3 

710 

58 Gly 

66 Tly 

8.8 x 10~ 4 

150Gy 

5.7 x 10 4 

< 1 

58 Gly 

2.3 Ply 

2.5 x 10“ 5 

200Gy 

2.0 x 10 6 

< 1 

58 Gly 

87Ply 

6.7 x 10~ 7 

250Gy 

7.4 x 10 7 

< 1 

58 Gly 

3.2 Ely 

1.8 x 10" 8 

500Gy 

5.4 x 10 15 

< 1 

58 Gly 

240 Yly 

2.4 x 10~ 13 

ITy 

3.0 x 10 31 

< 1 

58 Gly 

1.3 x 10 33 G% 

4.5 x IQ' 32 
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Table 3: Dependence on tr normalized to Riv^de) = Rev^de)- 


Turnaround 
(hr — hoc) 

a(t T ) 

Galaxies 
in introverse 

Riv(tr) 

Rev^t) 

f 

0 

0.75 

1.6 x 10 12 

39 Gly 

39 Gly 

1 

4Gy 

1 

10 12 

UGly 

52 Gly 

0.85 

10Gy 

1.6 

5.5 x 10 11 

49 Gly 

70 Gly 

0.70 

20Gy 

3.2 

8.8 x 10 1U 

53 Gly 

UIGly 

0.38 

30Gy 

6.6 

1.1 x 10 1U 

56 Gly 

290 Gly 

0.19 

40Gy 

13 

1.6 x 10 9 

57 Gly 

570 Gly 

0.10 

50Gy 

28 

1.4 x 10 8 

57 Gly 

1.2 Tly 

0.047 

lOOGy 

1.1 x 10 3 

2800 

58 Gly 

48 Tly 

1.2 x 10" 3 

150Gy 

4.3 x 10 4 

< 1 

58Gly 

1.9 Ply 

3.1 x 10~ 5 

200Gy 

1.5 x 10 6 

< 1 

58 Gly 

66 Ply 

8.8 x 10" 7 

250Gy 

5.5 x 10 Y 

< 1 

58 Gly 

2AEly 

2.4 x 10" 8 

500Gy 

4.0 x 10 15 

< 1 

58 Gly 

180 Yly 

3.2 x 10“ 16 

ITy 

2.2 x 10 31 

< 1 

58 Gly 

9.7 x 10 ' i2 Gly 

6.1 x IQ’ 32 


4 Bounce from Contraction to Expansion 

The CBE contracting universe contains no matter. It does, however, contain a known 
amount of radiation and in approaching the bounce this radiation dominates the energy 
and entropy. During the expansion era, at t = t 0 its energy contribution is D 7 (f 0 ) = 
1.3 x 10 -4 . Its entropy will be discussed below. 

The normal Friedmann equation without cyclic entropy is, from Eq.(j5]) and using the 
radiation equation of state 3 p = p 


a(t) 


8i tG 


p 7 (t), 


a(t ) 3 

and we need to calculate aits) which appears in this contraction version of Eq.( 

I^Tor(ts) — 1| = 


( 12 ) 


a(t 


B 


(13) 


Piif) 


Pji^o) 

a(t) 4 


Using 


(14) 

































in Eq. (lT2]) leads to the following 


da{t) _ 8vrGp 7 (t 0 ) m _ 3 
dt ~ 3 a[t) 


During the radiation-dominated era we have, from Table 1 that 



a(t) = C 7 f 1,/2 , 

(16) 

where 

C 7 = 2.3 x 10 _1 V 1/2 . 

(17) 

Now we can use 

O u 1 87rGp 7 (t 0 ) 

3 H(t 0 y 

(18) 

and 

\k\ = H{t 0 ) 2 

(19) 

to rewrite Eq.fllSl) 

in the more useful form 



S d{a) - S dt f3/2 

(20) 

and inserting 

\k\ = H(t 0 ) 2 

(21) 

in Eq.(TT3l) and performing the integrals gives the flatness result 



\^TOTAL{t b) ““ 1| — O^ts, 

(22) 


in which the coefficient, Cq, is readily calculated to be 


Ca = 




AH(t 0 ) 2 n^t 0 ) 


= 3.9 x 10 


-16 


-1 


(23) 


Eqs. (l22jl and (f23ll now give, ignoring prefactors which are 0(1), at the electroweak and 
Planck times, tEw = 10” 10 s and tpi anc k = 10~ 44 s, respectively 


| ^TOTAL (j^EW ) — 1 | ~ 0(10 26 ) 


(24) 


and 

| ^TOTAL(tpianck) ~ 1| ~ O(10 _6 °) (25) 

which are nothing more than Eqs. ([7]) and ([8]) respectively, now rederived using time 
reversal acting on the expansion era. 
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But the requirement of cyclic entropy importantly adds one more multiplicative factor 
on the right-hand-side of Eq. (l22j) because the CBE contracting universe is described by a 
modified scale factor a(t) = fa(t ), where f is the usually small fraction given in the last 
column of Tables 2 and 3 and in the fourth column of Table 4. 

With cyclic entropy, for the contracting universe, Eq. fllSl) must therefore be replaced by 

daif) = 87rGp 7 (f 0 ) _ 3 

dt 3 

Inserting now the defining CBE relationship d(t) = fa(t) we easily find that Eq. (l22|) 
inherits a multiplicative factor which is a fourth power of the CBE contraction fraction, 
/ 4 , giving the new flatness result 

|Ototal(Lb) — 1| = (27) 

Let us illustrate the remarkable implications of Eq. (l27)l by two examples, using the results 
displayed in Tables 2 and 3. 

As a first example, consider the turnaround time tr satisfying (tr — Ide) = 150Gly. This is 
the shortest time after which the visible universe contains < 1 galaxies and so is where the 
CBE assumption first becomes straightforward to implement. For this case / = 3.1 x 10~ 5 
and hence, ignoring prefactors which are 0(1), / 4 ~ 1CT 18 and therefore 


l^TOTAL^tEw) ~ 1 

- O(10~ 44 ) 

(28) 

and 

l^TOTAL^tPlanck) ~ 1 

~ O(10" 78 ). 

(29) 

As a second example, consider the turnaround time satisfying ( tr-toE) = 1 Ty so 
from Table 3, / 4 ~ 10~ 126 . This leads naturally to the remarkable initial conditions 

that, 

Ototal(Lew) — 1 

~ o(io- 152 ) 

(30) 

and 

\^TOTAL{tPlanck) ~ 1 

- O(10~ 186 ). 

(31) 

The total energy density at the present time t = 
density in these two examples. It is respectively 

to becomes extremely close to the critical 

Ototal^o) _ 1 

- O(10~ 18 ), 

(32) 
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for the choice [tx — t be) = 150Gh/,and 


I^total(O) — 1| ~ 0(10 126 ), (33) 

for the case (tx — t DE ) = 1 Ty. 

From a practical and observational point of view, however, the density perturbations of 
the CMB render it impossible to check Eqs. (l32p or (1331) to more than five decimal places. 


Homogeneity at the bounce is ensured by low entropy. Let us examine the value of the 
dimensionless entropy SxoTALits)/k which originates only from electromagnetic radiation, 
Stotal^b) = 

First note that at the present time the CMB in the visible universe has dimensionless 
entropy given in Weinberg’s book as |23| 

S,(t 0 )/k = g*V vu T 3 ~ 10 88 . (34) 


The radiation entropy remains constant during adiabatic contraction so the bounce entropy 
can be reliably estimated by using the extensive property that entropy is proportional to 
volume: 

S^(t B )/k ~ 10 88 /(f T ) 3 , (35) 

where we have displayed explicitly the important dependence of / on the turnaround time 
tr- 

Using Table 4, the bounce entropy for (tx — toE) = 150 Gy is thus S 1 (t B )/k ~ 10 74 which, 
while a very large number, is nevertheless extremely small compared to the present total 
entropy of the visible universe, Stotal^o) ~ 10 124 , and may be sufficiently small to ensure 
the observed homogeneity and isotropy. 

On the other hand, if S^(tB)/k ~ 10 74 is still insufficiently low, then as can be confirmed 
from the final column of Table 4, one may increase to a trillion years the turnaround time, 
{t T — t D E) = ITy, and thus render the bounce entropy identically zero, S^(tB)/k = 0, 
whereupon homogeneity and isotropy can become exact. 
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Table 4: Dependence of bounce entropy S 1 (ts)/k on (tr — 


Turnaround 
(tr — t DE ) 

a(t T ) 

Galaxies 
in introverse 

f 

Bounce 

entropy (S 7 /k) 

0 

0.75 

1.6 x 10 12 

1 

- 

4Gy 

1 

10 12 

0.85 

- 

10Gy 

1.6 

5.5 x 10 11 

0.70 

- 

20Gy 

3.2 

8.8 x 10 1U 

0.38 

- 

30Gy 

6.6 

1.1 x 10 1U 

0.19 

- 

40Gy 

13 

1.6 x 10 9 

0.10 

- 

50Gy 

28 

1.4 x 10 8 

0.047 

- 

lOOGy 

1.1 x 10 3 

2800 

1.2 x 10" 3 

- 

150Gy 

4.3 x 10 4 

< 1 

3.1 x 10~ 5 

~ 10 74 

200Gy 

1.5 x 10 e 

< 1 

8.8 x 10~ 7 

^lo 70 

250Gy 

5.5 x 10 7 

< 1 

2.4 x 10~ 8 

- 10 65 

500Gy 

4.0 x 10 15 

< 1 

3.2 x 10~ 16 

- 10 5U 

ITy 

2.2 x 10 31 

< 1 

6.1 x 10~ 32 

~ 0 


5 Discussion 

Previous work on cyclic cosmology, especially back in the twentieth century, was stymied by 
ignorance of the accelerated expansion which was first discovered only in 1998. The no-go 
theorem or Tolman Entropy Conundrum (TEC) implicitly assumed decelerated expansion 
and hence no extroverse surrounding the introverse into which to jettison the huge entropy 
which accumulates from irreversible processes during expansion. 

One prediction of the cyclic entropy model is that the presently observed flatness which 
can be expressed as |O(f 0 ) — 1| — 0 is precise, possibly being accurate to 18 or even 128 
decimal places depending on the turnaround time as exhibited in our Eqs. ()32l) and ()33lh 
This precision arises from the Come Back Empty (CBE) assumption (9j where the visible 
universe after turnaround generally has a radius only a very small fraction / of the radius 
of the previously expanding visible universe. Such theoretical precision far exceeds that 
of any possible observational measurement of O(to)- 

In summary, cyclic entropy provides an alternative to inflationary cosmology, at least as 
far as providing an alternative explanation of the observed homogeneity and flatness. To 
be a more complete alternative would require an explanation of the density perturbations 
which seed structure formation. In inflationary cosmology these are explained by quantum 
fluctuations of an inflaton field. In the alternative of cyclic entropy, such fluctuations might 
reasonably be expected to arise due to quantum effects during the contraction era. 
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